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Abstract 

The many-body Schrodinger dynamics of a one-dimensional bosonic Josephson junction is in- 
vestigated for up to ten thousand bosons and long times. The initial states are fully condensed 
and the interaction strength is weak. We report on a universal fragmentation dynamics on the 



to/) 



many-body level: systems consisting of different numbers of particles fragment to the same value 
at constant mean-field interaction strength. Using the Bose-Hubbard model we show how the value 
of the universal fragmentation dynamics can be predicted from the initial state. The analysis al- 
lows the prediction of the extent to which many-body effects become important in the dynamics at 
much later times. Even for the largest particle numbers and the weakest interaction strengths the 
dynamics is many-body in nature and the fragmentation universal. Implications are discussed. 
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Over the last few years ultracold bosons in double-well potentials have made it possible 
to explore quantum many-body dynamics in a highly controllable manner. Among the phe- 



nomena predicted or observed are tunneling and sel 



m , Elapse and revival seo.uer.ees fl, squeezing H matter wave interferon^ fl, frag- 
mentation and the approach to an equilibrium state |8Hl2|. While tunneling, self-trapping 
and Josephson oscillations have explanations on the mean- field level, other phenomena, e.g., 
collapse and revival sequences require many-body treatments such as the Bose-Hubbard 
model |l| or even the full many-body Schrodinger equation 



trapping |l-|5|, Josephson oscillations 



a, 
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In the context of bosonic Josephson junctions self-trapping - the inhibition of tunnel- 
ing due to interparticle interactions - is a theoretically and experimentally well-studied 
phenomenon jl-4|. Within two- mode Gross-Pitaevskii (GP) mean- field theory an explicit 
criterion was derived that determines whether self-trapping occurs for a given condensed 
initial state or not jl, 3|. The mean-field theory reduces the many-body quantum dynamics 
to that of a classical nonrigid pendulum and self-trapping occurs from a critical value of the 
interaction strength onwards (jj, [2 L The theory was successfully applied in the description of 
experiments at short time scales |3Hl|. However, for longer times the dynamics can leave the 
realm of mean-field theory and make a many-body description necessary, especially when 
the number of bosons is in the hundreds and the interaction strength is strong enough for 
self-trapping to occur, see, e.g., Refs. jl, Q, 14, 15|. 



Here, we investigate the long-time many-body dynamics of Bose-Einstein condensates 
(BECs) forming bosonic Josephson junctions with up to ten thousand bosons in the conden- 
sate. We consider interactions that are so weak that self-trapping cannot occur. The BECs 
are initially fully condensed. We then solve the time- dependent many-body Schrodinger 
equation for these BECs numerically. This letter contains the following main results. Firstly, 
we report on the existence of a universal many-body fragmentation dynamics in bosonic 
Josephson junctions. Fragmentation is a basic many-boson phenomenon which is - usually 
- strongly dependent on the number of bosons in the system js-12|. Here we find that 
systems consisting of different numbers of bosons all fragment to the same value at fixed 
mean-field interaction strength. The appearance of such universal fragmentation dynamics 
on the many-body level is unexpected. Moreover, this dynamical fragmentation implies the 
breakdown of GP mean-field theory. Above the critical interaction strength for self-trapping 
mean-field theory is known to break down on a time scale that increases logarithmically 



with the number of particles 



15] . Here, we find such a logarithmic breakdown even below 



the critical interaction strength. 

Secondly, using the Bose-Hubbard (BH) model we show how the universal fragmentation 
dynamics is related to the initial conditions. For a whole class of condensed initial states 
our analysis makes it possible to predict the fragmentation of the BEC after the collapse of 
the density oscillations, and thus to predict the extent to which the many-body dynamics 
for these condensed initial states remains within the realm of mean-field theory. Thirdly, we 
show that even for extremely weak interactions the dynamics remains many-body in nature, 
indicating that there is no limit at which GP theory becomes valid at long times. The results 
allow the extrapolation to larger particle numbers. 

In order to compute the time-evolution of the many-body Schrodinger equation we use 
the MultiConfi gur ational Time-Dependent Hartree for Bosons (MCTDHB) method see 
also Refs. 1 2 . T13I. \v\ fisj l . It is convenient to use dimensionless units defined by dividing 
the Hamiltonian by -Kp, where m is the mass of a boson and L is a length scale. The full 
many-body Hamiltonian then reads 

N 

H = Y,h(x l ) + Y,W(x i -x J ), (1) 

i=l i<j 

where h(x) = —\-^s + V(x) with a trapping potential V(x) and an interparticle interaction 
potential W(x-x') = \q5(x—x'). For the potential V(x) we choose a double-well constructed 
by connecting two harmonic potentials V±(x) = |(x ± 2) 2 with a natural cubic spline at 
\x\ = 0.5. Using the left- and right-localized orbitals 4>l,r constructed from the ground and 
first excited state of V(x) the parameters J = — {(f) L\h\(j) r) = 2.2334 x 10~ 2 , U = Xq J \4>l\ 4 , 
A = U(N - 1)/(2J), and t Rabi = n/J = 140.66 are computed. 

All one-body observables can be computed from the first-order reduced density matrix 
p^ l \x\x') = (^/\^^(x')^/(x)\^/} = Yli n[ af'(i)aP*(i')| where $f(x) denotes the bosonic 
field operator, the many-boson wavefunction, oq\x) the natural orbitals, and > 
rig 1 '', • • • the natural occupations with £\ n- 1 ^ = N. We suppress the time-argument whenever 
unambiguous. Bose-Einstein condensation is defined as follows: if an eig envalue np = 0{N) 
exists, the system is condensed 
said to be fragmented [8 
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Qf| . If there is more than one such eigenvalue, the BEC is 
2lfl . Here, it is convenient to define the fragmentation as 



/ = Y^ii>i n P /N. In GP theory / = 0% at all times, by construction. On the many-body 
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Schrodinger level this need not be the case and 100% > / > 0%. Fragmentation (/ > 0) is 



accompanied by a loss of coherence 



2l|. 



We compute the many-body dynamics for condensed states which are the GP ground 
states of the potential V + (x). Thus, initially the BECs are located in the left well. For 
each particle number N we choose the interaction strength A such that A = Xq(N — 1) 
is constant. The parameter A appears in the GP equation and is known as the mean-field 
interaction strength. Increasing N at constant A implies a decreasing interaction strength 
Ao- For the following computations we fix the mean-field interaction strength at A = 0.152 
which is equivalent to A = 1.33, i.e., well below the critical value A c = 2 for self-trapping 

13. 

Fig. [1] (top) shows the probability in the left well pz(t) = -h p^(x\x]t)dx as a 
function of time for N = 100 — 10000 bosons. The density tunnels back and forth through the 
potential barrier and eventually the density oscillations collapse. For larger N the collapse 
takes longer. The collapse occurs on the many-body level, but not within GP mean-field 
theory |l|, [l2|. We will therefore investigate the many-body nature of the BEC during the 
collapse in more detail. In Fig. [1] (bottom) the corresponding natural occupations are 
shown. Since the many-body wavefunction is initially condensed, there is only one natural 
occupation nf" — N at t — 0. However, as time increases a second natural orbital becomes 
occupied and the system starts to fragment. The nature of the BEC changes from condensed 
to fragmented as the density oscillations collapse. We stress here that these results represent 
the many-body Schrodinger dynamics of the Hamiltonian ([TJ) j^ . 

Looking at Fig. [T]we observe that the fragmentation dynamics is universal: at constant 
A systems consisting of different numbers of bosons fragment to the same value during the 
collapse of the density oscillations. At the end of the collapse the two largest natural occu- 
pations reach plateaus at about /N = 66% and /N = 34%. Thus, the fragmentation 
after the initial collapse of the density oscillations, f co i, settles at a value of f co i = 34%, 
regardless of the number of particles involved. Note that the GP mean-field has universal 
dynamics firmly built into the theory; the GP dynamics is identical for all systems with the 
same value of A. However, here we report on a universal fragmentation dynamics which is 
only present on the many-body level. Fragmentation is usually strongly dependent on the 



number of particles in the system 8Ml2|. Hence, the appearance of a universal fragmentation 
dynamics is unexpected. 



Let us discuss the time scales involved in this universal fragmentation dynamics. We 
define the fragmentation time Tf rag as the first time at which a certain fragmentation / is 
reached. Fig. [2] (top) shows the fragmentation time Tf rag as a function of TV. Clearly, Tf rag 
increases with N. For any value of the fragmentation, Tf rag is well described by a fit to the 
function T(N) = aln(l + bN), as is shown here for the values / = 4%, 9%, ...,34%. Thus, 
Tfrag grows logarithmically with N. Consequently, the fragmentation does not decrease or 
even disappear in the limit of large N. Even for N = 10000 bosons the fragmentation rises 
to about 10% in less than a dozen t Ra bi. 

Since GP theory does not allow BECs to fragment at all, Tf rag defines a measure for the 
breakdown of mean-field theory. It is well-known that mean-field theory can break down, 
especially for small particle numbers and on time scales of several tj? a f«, see, e.g., Refs. 
[iL Q; 14, 15 1 . In Ref. 15 1 it was shown that the time after which GP theory breaks down 
scales logarithmically with N for A = 4, way above the critical interaction strength for 
self-trapping A c = 2. Here, we find that the fragmentation time Tf rag scales logarithmically 
with N at A = 1.33, i.e., well below A c . Hence, the logarithmic breakdown of mean- field 
theory cannot be attributed to a mean-field instability |l5| and more research needs to be 
done to explain this scaling. 

So far we have established the universality of fragmentation of the many-body Schrodinger 
dynamics and discussed some implications. We will now analyze this dynamics on the basis 
of a simple model and show that the universality of fragmentation is a general phenomenon 
that exists for a wide range of initial conditions and interaction strengths. As a model we 
use the two-mode BH Hamiltonian 

Hbh = -J (b\b R + b R b L ) + | (b\V L b L b L + V R b R b R b R ) > (2) 

where the operators b^ L R create bosons in the orbitals 4>l,r, respectively. The probability in 
the left well reads pi(t) = jf(b\hi,)(t) and the fragmentation of the BEC reduces here to 
the occupation of the second natural orbital, / = /N < 50%. We consider the family of 
condensed states given by 

|tf > = 4= (VpW)b[ + VpW)b R ) N |0>, (3) 

where p R (0) = 1 - p L (0). 

We choose parameter values corresponding to A = 1.33 and solve the dynamics of the 
Hamiltonian (j2J). We find that for the entire family of initial states |\&o) with < pz,(0) < 1, 



BECs consisting of different numbers of particles fragment to the same value j co \ during the 
collapse of the density oscillations. Thus, also the BH fragmentation dynamics is universal. 
This can be seen in Fig. [3] (left) where the fragmentation / is shown as a function of time for 
N = 1000 and 10000 bosons. The fragmentation of initial states with pz,(0) = 1.0 reaches a 
plateau at about f co i = 32%, not far from the many-body Schrodinger result 34%, see Fig. [T] 
(bottom). Similarly, initial states with pi{0) = 0.8 fragment towards the value f co i = 16%. 
Thus, the universal fragmentation dynamics is a general phenomenon which exists for a 
whole family of initially-condensed systems. 

We will now point out a striking relation that allows the prediction of the value of f co i 
from the initial state. To this end we solve the BH eigenvalue problem Hbh\E u ) = E n \E n ). 
Fig. H2 (right) shows the fragmentation /„ of each eigenstate \E n ) as a function of E n /N. 
Note that the curves for the eigenstate fragmentation f n = f n (E n /N) practically do not 
depend on the particle number, as shown here explicitly for N = 1000 and 10000 bosons. 
Also shown are the expansion coefficients C n = (E n \^ ) of the initial states together with 
their energies per particle E/N = (^o\Hbh\^o)/N. As a function of E n /N the coefficients 
are approximately binomially distributed with a width ^iV" 1 / 2 . We find that for eigenstates 
with energy E n rs E the eigenstate fragmentation equals the fragmentation of the many- 
body system after the collapse of the density oscillations, namely f n (E n /N « E/N) = / co ;. 
This surprising empirical result links the condensed initial states (|3]) to the fragmentation 
of the many-body wavefunction after the collapse of the density oscillations, i.e., at much 
later times. We have checked this relation for up to N = 10000 bosons. However, the results 
shown in Fig. (J3]) suggest that an extrapolation to any number of bosons can be made. 
Based on these findings it is possible to predict to what extent the condensed initial states 
([3]) will fragment during the collapse of the density oscillations. We note that fragmentation 
manifests itself in observables such as the fringe visibility and corrrelation functions which 



are experimentally accessible, see, e.g., Refs. |l5l. l2lL l23|. 

Let us briefly touch upon the question whether fragmentation disappears for ever weaker 
interaction strength. Therefore, we choose an interaction strength which is only a tenth of 
the interaction strength considered so far, A = 0.0152 (A = 0.133). The results are shown in 
Fig. HI Again, the BECs fragment and the fragmentation dynamics is universal. Note that 
the time scale is much larger now. Initial states with Pl(0) = 1.0 now fragment to about 
fcoi — 48% which is higher than for stronger interaction, initial states with Px(0) = 0.8 
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to about f co i = 10.5% which is lower than for stronger interaction. These results are in 
agreement with the predictions of the here proposed model for universal fragmentation: we 
find again that f n (E n /N E/N) = f co i. The results also show that the GP mean- field 
dynamics is not recovered despite the weaker interaction strength. 

Summarizing, we have shown that there is a universal fragmentation dynamics in bosonic 
Josephson junctions at weak interaction strengths by solving the time-dependent many- 
boson Schrodinger equation. We have shown how the value of this universal fragmentation 
can be predicted from the initial state using the two-mode BH model. The GP mean-field 
dynamics is not recovered in the limit of large particle numbers, despite weak interaction 
strength. 
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10 t/t Rabi 20 30 

Figure 1. (color online): Universality of the many-body Schrodinger fragmentation dynamics. 
Shown is the probablility in the left potential well pz,(t) (top) and the corresponding natural occu- 
pations (bottom) as a function of time for different particle numbers. During the collapse of the 
density oscillations (top) the nature of the BEC changes from condensed to fragmented (bottom). 
The fragmentation reaches a plateau at / co ; = 34% for all particle numbers. Interaction strength: 
A = 0.152 (A = 1.33). Particle numbers: N=100 (blue), 1000 (magenta), 10000 (orange). The grey 
lines represent particle numbers in between, N=200, 500, 2000 and 5000. All quantities shown are 
dimensionless. 
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Figure 2. (color online): Fragmentation times. Shown are the times Tf rag that are needed to reach 
a given fragmentation / as a function of the number of bosons for the same parameters as in Fig. [TJ 
The data points (crosses) of all fragmentation values are well described by a logarithmic fit function 
(lines), here shown for the values / = 4%, 9%, . . . , 34%. This logarithmic dependence of Tf rag on 
the number of particles implies the breakdown of GP mean-field theory even in the limit of large 
particle numbers. Note that the logarithmic dependence occurs already at an interaction strength 
A = 1.33, well below the critical mean-field interaction strength A c = 2. See text for details. All 
quantities shown are dimensionless. 
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Figure 3. (color online): Universality of the BH fragmentation dynamics. Interaction strength 
A = 0.152 (A = 1.33). Left: Shown is the fragmentation dynamics of fully condensed initial states 
with different numbers of particles in the left well Pl(0). Parameter values: N = 1000, £>l(0) = 
1.0 :;N= 1000,p L (0) = 0.8 : 0, N = 10000,p L (0) = 1.0 : N = 10000,p L (0) = 0.8 : ■ The 
fragmentation reaches a plateau at about j co \ = 32% (/ co ; = 16%) for initial states withpx(O) = 1.0 
(pl(0) = 0.8). Right: Shown is the fragmentation f n of each eigenstate as a function of its energy 
per particle E n /N for different particle numbers (N=1000: 0, N=10000: □). Also shown are the 
energies per particle of the initial states E/N (vertical lines) together with the distributions of the 
expansion coefficients C n as shaded areas (N=1000: blue, N=10000: red). For eigenstates with 
E n ~ E the fragmentation of the many-boson wavefunction after the collapse is given by f co i = f n 
for any value of Pl{®)- All quantities shown are dimensionless. 
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Figure 4. (color online): Same as in Fig. [U but for ten times smaller interaction strength A = 0. 
(A = 0.133), far below the critical interaction strength A c = 2. See text for more details, 
quantities shown are dimensionless. 
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